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Abstract 

An oscillatory pattern in the smoothed quantum spectrum, which is unique for single- 
particle motions in a reflection-asymmetric superdeformed oscillator potential, is investi- 
gated by means of the semiclassical theory of shell structure. Clear correspondence be- 
tween the oscillating components of the smoothed level density and the classical periodic 
orbits is found. It is shown that an interference effect between two families of the short 
periodic orbits, called supershell effect, develops with increasing reflection-asymmetric 
deformations. Possible origins of this enhancement phenomena as well as quantum signa- 
tures of period-multipling bifurcations are discussed in connection with stabilities of the 
classical periodic orbits. 

§1. INTRODUCTION 



Possible occurrence of instability of superdeformed (SD) nuclei having the prolate shape with 
the axis ratio approximately 2:1 against the octupole-type reflection asymmetric deformation is 
one of the current topics of growing interest in high-spin nuclear structure physics. Regions in 
the {N,Z) plane where we can expect existence of reflection- asymmetric SD nuclei have been 
investigated mainly by means of the Strutinsky-type calculations of the collective potential 
energy surface (see also Refs. and |0] for other approaches). Concerning the physical condi- 
tion for the occurrence of the octupole instability, Nazarewicz and Dobaczewski S have recently 
discussed the connection between the dynamical symmetry of the anisotropic harmonic-oscillator 
with frequencies in rational ratio and the multi-cluster configurations. They have suggested 
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that the closed-shell configurations in the prolate SD oscillator potential, defined as having the 
frequency ratio Lo^jio^ = 2, might be unstable (stable) against the octupole-type refiection asym- 
metric shapes when the single-particle levels are filled up to the major shells with iVsh=even 
(odd), Ai'sh being the shell quantum number defined by Ng^ = 2n± + (see also the previous 
work, Ref. 0). Their suggestion is in good qualitative agreement with the realistic shell-structure 
energy calculation by Holler and Aberg. We call the A'sh-dependence of the octupole instability 
"odd-even effect in A^sh" • 

We have suggested in Refs. [jlO|, fill a possible relationship between the odd-even effect 
mentioned above and the "supershell effect" in refiection-asymmetric SD potentials. The general 
concept of supershell was originally introduced by Balian and Bloch@ in relation to the semi- 
classical theory of shell structure. Quite recently, the supershell effect has been observed, for 
the first time, in the mass abundance spectra of metal clusters. Theoretical analysis of this 
phenomena has been made by Nishioka, Hansen and Mottelson. 

As is well known, clustering of eigenvalues, that is, oscillating pattern in the energy-smoothed 
level density for single-particle motions in the mean field is called shell structure. ^ In the semi- 
classical theory, classical periodic orbits having relatively short periods are responsible for the 
clustering of the levels; the frequencies in the level density oscillation is determined by the cor- 
responding periods of classical motion. When two families of short period orbits interfere 
and produce an undulating pattern in the oscillating level density, this pattern is called super- 
shell structure, ii'ii In the case of the metal clusters, a beautiful beating pattern enveloping 
individual shell oscillations which is caused by the interference between the triangular and square 
orbits of an electron in a spherical Woods-Saxon potential has been demonstrated @'0 to nicely 
correspond to the experimental data. 

In the case of the SD nuclei under consideration, an interference effect is expected to arise 
between the classical periodic orbits with period T ^ 2tt/uj_[_ and those with T ^ 27r/ci;2 of a 
nucleon in the reflection-asymmetric SD potential. The main purpose of this paper is to show 
that the interference effect brings about another example of the supershell structure, which is 
intimately connected with the odd-even effect in A^sh mentioned above, and which is relevant to 
experimental investigations. It should be emphasized here that, contrary to the case of spherical 
potentials, the Hamiltonian describing the single-particle motions in a deformed mean fleld is 
non-integrable in general. Accordingly, our Hamiltonian system is a kind of mixed system whose 
phase space is composed both of regular and chaotic regions. As a consequence, properties of our 
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phase space change in a quite sensitive manner when the shape of the nuclear surface is varied. In 
this paper, we shall show that the supershell structure becomes more significant with increasing 
octupole deformation. Possible origins of this enhancement will be discussed in relation to the 
change in the properties of the classical periodic orbits as a function of the octupole deformation 
parameter. 

After briefly reviewing the semiclassical theory of shell structure in §2, we first apply in §3 
both the torus quantization method and the periodic-orbit quantization method to the case of 
the prolate SD oscillator potential. In this integrable limit, the supershell effect can be treated 
analytically. In §4, a reflection-asymmetric SD potential model is introduced and the supershell 
pattern in the quantum level spectrum calculated for this potential is exhibited. In §5, we inves- 
tigate properties of classical motions in this potential, like stabilities and bifurcation phenomena 
of the periodic orbits. In §6, we show that a nice correspondence holds between peak positions of 
the Fourier transform of quantum spectrum and periods of classical closed orbits; relative heights 
between peaks change as functions of the octupole- deformation parameter, providing us with a 
semiclassical interpretation of the origin and the octupole-deformation dependence of the super- 
shell structure. Here, quantum signature of the bifurcations is also discussed. Summary of this 
work is given in §7. 

A preliminary version of this work was previously reported in this journal. 

§2. SOME ELEMENTS OF SEMICLASSICAL THEORY OF SHELL STRUCTURE 

In this section we briefly review some basic elements of the semiclassical theory of shell struc- 
ture, which are necessary for later discussions. 

2.1. Torus Quantization 

To begin with, let us consider the case of multi-dimensional, integrable Hamiltonian system, 
where the Hamiltonian can be written as a function of only action variables Jj, being indepen- 
dent of angle variables 6i conjugate to them. Semiclassical quantization condition valid for such 
systems has been formulated by Einstein-Brillouin-Keller, and called torus quantization or the 
EBK quantization; 
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Ii{En,,..,nf) = f^pdq = n{n, + ai/4) , i = l,...,/ (2.1) 

where indices i represent mutually independent paths on /-dimensional torus constructed by clas- 
sical trajectories, are Maslov indices related to the singularities of the Van Vleck determinant 
appearing in the semiclassical propagator along the path i. Thus, the semiclassical level density 
is given by 



{"■} 

The summation on the r.h.s. may be rewritten using the Poisson sum formula into the form of 
topological sum over periodic orbits, il In Ref. |^5|, spherical systems are analyzed and clear 



correspondence between the topological sum and the periodic orbits is shown. In the spherical 
case, all periodic orbits satisfy the resonance condition, i.e., the frequency ratio of radial and 
angular motions are the same as that of topological indices. We shall apply this method to the 
SD harmonic oscillator potential in §3, and discuss the correspondence between the topological 
indices and periodic orbits. There, it will be shown that some "partially-resonant" terms play an 
important role giving rise to the supershell effect. 



2.2. Periodic-Orbit Quantization 

Next, let us consider the case of multi-dimensional non-integrable Hamiltonian system. For 
such systems, as is well known, the periodic-orbit quantization method provides us with a useful 
base toward understanding the correspondence between classical periodic orbits and properties 
of quantum spectra. ©''B This theory is essentially based on the path integral formalism of 
quantum mechanics. The first step is to express the quantum level density g{E) = J2n ^{E — En) 
as a trace of the energy-dependent Green function; 

1 1 
g{E) = ImTr- 



vr E + ie-H 



--Im / dqG{q,q;E) . (2.3) 

71 J 



The Green function G{q" , q'; E) is a Fourier transform of the transition amplitude K{q" , t; q', 0) = 
{q"\exp{—itH/h) \q')-9{t), and we can express it in the path integral form. Evaluation of the 
path integral by the stationary phase approximation (SPA) extracts the classical trajectories. 
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The Fourier transformation is also performed by means of the SPA. Finally, the trace integral 
appearing in Eq. ( p.3|) extracts the periodic orbits and one obtains the following expression called 
the Gutzwiller trace formula: 



g{E) ^ g{E) + ^ A„^(E) cos {nS^{E)/h - (7r/2)/i„,,) , (2.4) 



n,7 



where g{E) denotes the average level density and the second term on the r.h.s. represents the 
oscillating part. The summation is taken over all periodic orbits and their multiple traversals. 5*^ 
is a classical action along the orbit 7, = §^p-dq, and fij is a Maslov phase. The amplitude 
factor Any depends on the phase space structure about the periodic orbit 7, as we shall discuss 
in §2.4. For sufficiently isolated orbits, the trace integral is well approximated by the SPA and 
the amplitude factor for the n-fold traversal of orbit 7 can be written as© 

1 

"^"^ ~ ^ J|det(l -i^ ' ^^'^^ 



where and represent the period and the monodromy matrix of the primitive orbit 7, 
respectively. This expression is known to work well for chaotic systems such as billiards. 



2.3. Stability of Classical Trajectories 

The amplitude factor in the trace formula is related with the properties of phase space around 
the periodic orbits. Let us write the Hamilton equation in 2/-dimensional phase space as 

= AVH , (2.6) 

dt ^ ^ 

with 




and consider the time evolution of the deviation 5Z{t) from the reference classical trajectory 
Zo(t). To the first order in SZ, we obtain 

^5Z = m5Z , (2.7) 
dt ^ ' 

where Ti is the Hessian matrix defined by 
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Knowing Hit), Eq. can be easily integrated; 



= exp 



A f &t' n{Zo{t')) 
Jo 



SZ{0)=Mit)SZ{0) . (2.9) 



Ai is called the stability matrix. It is real and symplectic; AA^-^A^^ = A4^^. When one take 
a periodic orbit as the reference trajectory and the period T as time t, the stability matrix is 
particularly called monodromy matrix = A4{T^).^ It is known that the eigenvalues of 
are independent of the initial condition Z{0) on the orbit. According to the symplectic property 
of the monodromy matrix, its eigenvalues appear in pairs (+/— )(e°', e~°), where a is real or pure 
imaginary. When a is pure imaginary {a = iw), the orbit is stable and torus exists surrounding 
it. is called winding number of the torus. When a is real {a = AT), the orbit is unstable and 
A is called the Lyapunov exponent which measures the degree of instability. 



2.4. Magnitudes of Shell Effects 

Strength of the shell effect depends mainly on three factors, to be discussed below, associated 
with the periodic orbits. El 

The first factor is the degeneracy of the orbit. Here the term 'degeneracy' means the number 
of independent continuous parameters (additional to energy) that specify a certain orbit from 
continuous family of orbits having the same action. For example, planar orbits in a spherical 
potential form continuous family generated by rotation and the degeneracy is generally 3, since 
a certain orbit belonging to this family is specified by three Euler angles. As illustrated by the 
above example, degeneracy is related to continuous symmetry of the system. These degeneracies 
correspond to the unit eigenvalues of M. 

The second factor is the stability of the orbit. For non-integrable systems, evaluating the trace 
integral by the SPA, one sees that the amplitude behaves as 

oc , ^ , = , (2.10) 
V'|det(l-M«)| 

where M is a reduced monodromy matrix in which degrees of freedom corresponding to the unit 
eigenvalues of M are excluded out. The more unstable is the orbit, the weaker is its contribution 
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to the level density, because it has a large Lyapunov exponent and the denominator on the 
r.h.s. becomes large. The above proportionality is valid only when all the eigenvalues of M are 
sufficiently distant from unity. However, one of the eigenvalues may happen to be very close to 
unity. This is called nonlinear resonance where two frequencies of the torus coincide with each 
other and gives rise to a periodic orbit bifurcation. Namely, the period n-upling bifurcation occurs 
when det(l — M") = 1. In this resonance region, one have to use a more sophisticated treatment 
than the SPA; for example, the uniform approximation using the resonant normal form. Such 
a procedure is formally discussed by Ozorio de Almeida and Hannay, il but, to the best of our 
knowledge, application of this theory to mult i- dimensional, non-integrable Hamiltonian system 
has not been performed yet. 

The third factor is the phase space volume occupied by an orbit. It is not important in our 
analysis because it is insensitive to the variation of potential parameters. 

Let us examine how these three factors enter in the amplitude factors for different types of 
periodic orbits. First, consider a chaotic orbit, that is, a well- isolated orbit whose degeneracy 



equals zero. Its amplitude factor is given by Eq. (2^). There appears the same stability factor as 
Eq. (p.lO|) , and the period T measures the phase space volume of the orbit. Next, as an example 
of non-isolated orbits, let us consider orbits in axially symmetric deformed potentials, whose 
degeneracy equals one corresponding to the rotation about the symmetry axis. It is convenient 
to use the cylindrical coordinates (p, <y9, z). One should then perform the integral uniformly with 
respect to in the trace, because in this direction periodic orbits exist continuously. Thus we 
obtain the following expression of the amplitude factor for these orbits (see Appendix A): 

^n.{E) = , , (2.11) 



B^= I dt 



7 
-1/2 



Here, contains the ffist and the third factors mentioned above. The period n-upling bifurcation 



occurs when TrM" = 2. 



§3. SUPERSHELL EFFECT IN THE SD OSCILLATOR 

In this section we apply the semiclassical theories to the axially-symmetric 2:1 deformed 
harmonic oscillator Hamiltonian 
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j:^-m, (3.1) 

i=x,y,z 

where ojx = ojy = = 2uJz = 2uJsh, and we discuss how the supershell structure emerges in 
this case. We compare the two semiclassical quantization methods summarized in the preceding 
section, and discuss their relations. 



3.1. The Periodic Orbit Method 



In this subsection we analyze the supershell effect in the SD oscillator defined by (|3.1|) using 
the Gutzwiller trace formula. The trace formula can be derived also for such an integrable system 
if the degeneracy (mentioned below) is properly taken into account. According to Ref. the 
semiclassical level density may be written as 



The first term on the r.h.s. represents the mean level density; 



(3.2) 



K^^sh) 



(3.3) 



The second and the third terms are the oscillating parts representing the shell effects. The 
superscript (II) and (IV) denote the degeneracies of the periodic orbits, g^^^^'' (E) is a contribution 
from four-hold degenerate orbits whose periods are multiples of T'^^' = 27r/co'sh. The present 
model is very special in the sense that all trajectories are periodic, and one should explicitly 
perform four integrals (corresponding to the degeneracy) in evaluating the trace formula. Thus 
one obtains the following expression; 

E^ 



8(^u;sh)- 



cos 



oo 



^(IV) 

m I — (4 + 4 



7r 



E 



m=— oo 



■ cos 



m 



N 



h 

/^7-(iv) 
htOsh ( N 



5 
2 



9{E) . 



(3.4) 
(3.5) 



where S**^^' = ET^^^^ is the action integral along the primitive periodic orbit and the sum over 
m accounts for multiple traversals. The last expression is obtained by using the Poisson sum 
formula 
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exp(27rimA) = ^ 6{A-N) . (3.6) 

m=— oo N=—oo 

In Eq. ( p.2| ), g^^l{E) is the contribution from two- fold degenerate orbits whose periods are odd 
integer times T*"' = It is obtained in a similar manner, except that the integrations with 

respect to the z direction may be performed by the SPA. The result is written as 



2^ — nrCOS 



r7i=— oo 



K^^sh)^ 



(2m + i; 



(3.7) 



E(-)"^^^4^-^-shfiV+^)) , (3.8 



N 



2. 

when the sum over m' accounts for multiple traversals of the primitive periodic orbit. The 



expressions ( |3.4|) and ( pTTl) were first derived in Ref. |29 



Summing up the above three contributions, we obtain the degeneracy of the A^-th shell as 

([Ar/2] + l)([Ar/2]+2) 3 
dN = + - , (3.9) 

where [*] is the Gauss symbol. The first term on the r.h.s. correspond to the exact degeneracy of 
the quantum spectrum. We thus see that the result obtained by the trace formula is very accurate 
(the deviation from the exact quantum result is only 3/32). 

Now let us focus our attention on a smoothed density of levels with finite energy resolution 
6E = huJsh- It is then sufficient to consider finite number of periodic orbits of short periods 
satisfying the following uncertainty relation; 

27rh 

T < T^ax = ^ . (3.10) 
oE 

As far as gross properties of the level density is concerned, therefore, the well-known problem of 
the long time propagation in the semiclassical approximation does not occur. 

We show in Figs. |l]-a) and |l]-b) the contributions from the families of periodic orbits with 
periods 27r/co'sh and 27!'/uj±, respectively, and in Fig. |l|-c) the sum of them. There appears an 
undulating pattern in the level density due to the interference of the above two families of periodic 
orbits, which is just the supershell structure. Thus, one sees that the supershell structure emerges 
from this interference effect. 
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3.2. The EBK Method 



Defining the action-angle variables (Jj, 6i) by 



Pi = \ 2 Muili sin 9i , 

^' _^ (3.11) 



Qi = ^2Ii/Mui COS 9i , 

the Hamiltonian ( p.lD is written as a function of only action variables as Hq{I) = uj-I. The 
Maslov indices are 2 for all paths and therefore the EBK quantization condition becomes 

h = him + 1/2) . (3.12) 

In the present case, this EBK quantization gives exact quantum eigenvalues: E = J2i huJi{ni+l/2). 
Now let us investigate the roles of classical periodic orbits in giving rise to the supershell structure 
in the quantum spectrum. For this purpose, we use the method of topological sum. il The 
semiclassical level density is written as 

g{E) = Y,6[e -J2huJi{ni + a,/A)) . (3.13) 



Using the Poisson sum formula ( |3.6|) , one can rewrite Eq. ( |3.13 ) as 



g{E) = go{E)+J2g,,{E), (3.14) 

where 

gAE) = ^ rdISiE-Lj-I)expi2mM-iI/h-(x/A)) , (3.15) 
n Jo 

and the summation is taken over all the combinations of integers, M = {M^, My, M^). Here, 
go{E) represents a mean level density corresponding to the Thomas- Fermi approximation; 

9.(£) = J,f d/^(£-„./) = ^. (3,16) 

The remaining terms with non-zero M represent the oscillating part responsible for the shell 
structure. To simplify the expression, we introduce the notation fi = 2'KMiUJsh/^i- The dominant 
contribution comes from terms satisfying the resonance condition M = M* oc lo, i.e., fx = fy = 
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fz (in the present case, M = m(2,2,l)). Carrying out the integration with respect to / and 
denoting the sum over such resonant terms as gj^^^, we obtain 

9liE) = — ^ Y: e^— (^Z'^---^/^) . (3.17) 

Next, let us consider the 'partially resonant' terms which satisfy the condition fi = fj = fj_ ff.. 
Carrying out the integration with respect to I, they are evaluated as 

^A/ (part, res.) 



4(^^sh)2 w^-fk) if±-fkr 

1 E _^if_LE/hui^i,^inM-a/2 ^ (3.18) 



where the second term on the r.h.s. is neglected because it is higher-order in h. Let us then take 
terms with fx = fy = f±^ fz, and write M as (m', m', I). Summing over terms with odd-m' and 
arbitrary integer /, we obtain 

1 E 



2TTim' E /2hiOsh , ^-2TTi{2m' +1} /2 



m'=odd l=~oo 



rp CO 

Y 6*''^^'"+^^^^/^'^='^"^/^^ . (3.19) 



8(^a;sh) 



2 

m=— oo 



The last expression is obtained using the expansion formula of cosecond in partial fractions: 

oo ^ y 

z — Iti 



cosec^ = ■ ^^-^O) 

(=— oo 



It can be easily shown in a similar way that the sum of other terms which are the same order in 
h as ( p.l9|) vanishes. 



3.3. Relation between the Two Methods 

Now, let us discuss the correspondence between the results obtained by the periodic-orbit 
method and the EBK method. It is evident that the two terms g^^^^{E) and g^^^l{E) in the 
EBK treatment are identical with the contributions g^^J{E) and g^^l{E) which are evaluated 
by the trace formula for the periodic orbits with periods 2Tr/uJsh and 27:/u±, respectively. This 
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result is very instructive to understand the physical meanings of the resummation with respect 
to the indices M by the use of the Poisson sum formula. The indices satisfying the resonance 
condition, M = m(2, 2, 1), correspond to a family of classical orbits with periods 2m7r/c(Jsh. One 
can examine this by comparing Eqs. ( p.4|) and (|3.17|) . On the other hand, for partially- resonant 
contributions, we find the correspondence in the following way. Comparing Eqs. ( |3.7| ) and ( ^.19[ ), 
we notice that the family of planar orbits in the {x, y) plane corresponds to the summation over 
indices M = (1, 1,/) with — oo < / < oo. Its m-fold traversals are related with M = {m,m,l). 
These partially-resonant terms play important roles in formation of the supershell structure in 
the present model. 



§4. REFLECTION-ASYMMETRIC SD OSCILLATOR MODEL 

4.1. Model Hamiltonian 

Let us consider a model Hamiltonian consisting of an axially-symmetric 2:1 deformed harmonic 
oscillator and a doubly-stretched octupole (I30) deformed potential; 

2 / 2 \ " 

H = ^ + Mu;l^-^-\,or'Y,oi9)j , (4.1) 

where double primes indicate that the variables in parenthesis are defined in terms of the doubly- 
stretched coordinates q'/ = (ui/ujQ)qi, and = (cUxCOyCUzY^^ . As emphasized by Sakamoto 
and KishimotoB, the doubly-stretched coordinates are suited to description of systems having 
quadrupole equihbrium deformations, and possess several advantages over the usual coordinates; 
for example, the center of mass motion is exactly decoupled from the octupole-type deformations 
described by the above Hamiltonian. Note that the doubly-stretched octupole operator is in fact a 
linear combination of the ordinary dipole and octupole operators, although we sometimes omit the 
adjective "doubly-stretched" for simplicity. In (|4.1| ), we adopt the quadratic radial dependence 
for the octupole-deformed potential for the ease of taking into account the volume conservation 
condition. By requiring the volume surrounded by an equipotential surface to be independent of 
the octupole deformation parameter A30, the Aao-dependence of Uq is determined as 



'^o(A3o) = ^o(O) 



1 



-3/2 



1/3 



(4.2) 
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We note that the average level density g{E) is independent of A30 when uq satisfies Eq. ( [4 .21) . Let 
us define dimensionless variables as 



Pi ^MhujQ Pi , 

^ ^h/MiUo q, , (4.3) 
H ^ huoH . 



Then the Hamiltonian ( [4.1|) becomes 

2/2 \ " 

H=^+[j-hor'Yso{9)j . (4.4) 

Since this potential is a homogeneous function of the second order in coordinates, the scaling 
relation 

H{ap, aq) = a^H{p, q) (4.5) 

holds. Thus, if {p{t),q{t)) is a solution of the Hamilton equation with energy E, {ap(t) , aq{t)) 
is also a solution but with energy a'^E. Namely, once the classical properties of the system are 
known on a certain energy surface Eq, properties on other energy surface E are obtained by scale 



transforming the phase space variables {p,q) to {ap,aq) with a=jE/EQ 



4.2. Supershell Structure 



Figure |^ shows the oscillating part of the level density for the Hamiltonian ( |4.4| ) calculated by 
means of the Strutinsky method. A characteristic property of the oscillating level density is that 
it exhibit the supershell pattern. Figure ^ gives a phenomenological illustration of the concept 
of the supershell. It is seen from this figure that the oscillating level density can be represented 
as a superposition of trigonometrical functions, cos{ET^/h) with ^ 2tt/uj±_ and 27r/c(Jsh, in a 
good approximation. We shall confirm in §6 that this supershell pattern in fact arises from an 
interference effect between two families of classical periodic orbits with periods ~ 27: /us^ and 
27r/c(jj_. 

It should be recalled here that the important factor from the point of view of gaining the shell- 
structure energy is not the heights of the maxima but the depths of the minima in the oscillating 
level density. Needless to say, the minima in Fig. |^ correspond to the closed-shell configurations 
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with respect to the SD major shell quantum number Nsh- Let us notice how the depths of the 
minima change as functions of A30. Then we find that the minima associated with the odd-A^^sh 
closures becomes shallower as A30 increases, whereas those at the even-A^'sh closures are tough. 
Consequently, the odd-even staggering of the minima with respect to the A^sh quantum number 
develops with increasing A30. Possible mechanisms of the enhancement of the supershell structure 
will be discussed in 56. 



§5. CLASSICAL ANALYSIS 

In this section, we discuss the classical-mechanical properties of the single-particle motion in 
the reflection-asymmetric SD potential defined in the preceding section. 



5.1. Poincare Map 

Let us examine classical phase space structure by plotting the Poincare map. Since our 
Hamiltonian is axially symmetric, it reduces to a two-dimensional one with the cylindrical coor- 
dinates (p, z) and with a definite angular momentum = m; 

H =^{pI + pI) +V,s{p,z-m) , (5.1) 

where 

V;.(.,p;P.) = ^ + - Xso^^ ^^=^ . (5.2) 

We can examine the Poincare map for each value of m. It is convenient to choose the Poincare 
section E as the surface with pp = 0, which is intersected by any trajectory. Figure § shows 
calculated Poincare maps {z,Pz) for the Hamiltonian ( [4.1|) with various values of the octupole 
deformation parameter A30. We see that the system is quasi- integrable for small A30 and almost 
all the phase space is foliated by KAM tori. With increasing A30, however, chaotic regions begin 
to spread out from the hyperbolic points. Figure ^ shows Poincare maps for different values of 
Py,. We note that the phase space volume corresponding to the (p, z) degrees of freedom contract 
and the system becomes more regular as p^p increases. 

Figure ^ shows Poincare maps for the surface of section {p,Pp) with z = 0. The origin 
corresponds to the linear orbit along the z-axis and the structure around it is exhibited. 
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5.2. Periodic Orbits and Their Bifurcations 



While all trajectories are periodic when A30 = 0, only very limited trajectories remain as 
periodic orbits when A30 7^ 0. In the Poincare section, centers of tori correspond to stable periodic 
orbits, while saddles to unstable ones. We calculate the periodic orbits by the monodromy method 
proposed by Baranger, Davies and Mahoney. In this method, periodic orbits are found in an 
iterative manner starting with approximate closed curves. By gradually changing A30, we can 
use the periodic orbits found for a slightly smaller value of A30 as inputs for this procedure. 
Figure shows short periodic orbits for Hamiltonian ([4.1|) with A30 = 0.4 obtained in this way. 
Also shown in Fig. are planar orbits for A30 = 0.3 ~ 0.4. As A30 increases, the phase space 
structure becomes more complicated due to bifurcations of stable periodic orbits. For example, 
a period-tripling bifurcation of orbit A occurs at A30 — 0.36. Thereafter, orbit A bifurcates into 
orbits 3A (triple traversal of orbit A), E and F. In the Poincare map (see Fig. ^) for A30 = 0.37, 
we can see three resonant island chains surrounding the central KAM torus, which are associated 
with the newly-born periodic orbits E and F. Likewise, a period-doubling bifurcation of orbit B 
occurs at A30 — 0.4, from where orbit B bifurcates into orbits 2B (double traversal of orbit B) 
and K. Many higher order bifurcations occur almost everywhere in regular regions of the phase 
space. Properties of the calculated periodic orbits are summarized in Table I for several values of 

-^30- 



§6. SEMICLASSICAL ANALYSIS 

6.1. The Cause of the Enhancement of the Supershell Effect 

As mentioned in §2.4, strength of shell effect associated with a periodic orbit is mainly de- 
termined by degeneracy and stability of the orbit. Let us discuss how these properties change 
when the octupole deformation is added to the SD oscillator. When A30 = 0, orbits with the 
period 2ti /uj^h and those with the period Itx joj^ have different degeneracies, 4 and 2, respectively. 
Therefore, the shell effect originating from the former families of orbits is much stronger than 
that from the latter. As a result, the interference effect between the two families of periodic orbit, 
i.e., the supershell effect, is rather weak. When A30 7^ 0, in general, degeneracies of the orbits 
reduce to 1. (For orbits D and A' in Fig. having special symmetry, the degeneracy is 0.) Thus, 
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generally speaking, shell effects at A30 7^ are expected to become weaker in comparison with 
those in the A30 = limit, and it seems hard to understand the enhancement mechanism of the 
supershell effect at A30 7^ in terms of the degeneracy property. 

Next, let us consider the other factor, i.e., stability of periodic orbits. In Fig. ^ we show 
calculated values of TrM for relevant orbits as functions of A30. At A30 = 0, orbits with the 
period 27r/c(Jsh are resonant (TrM = 2), while orbits with the period 27t/uj± are non-resonant 
and take TrM = —2. With increasing A30, TrM for orbits 2A(double traversal of orbit A) and 
B decrease and deviate from 2. At A30 — 0.4, a period-doubling bifurcation of orbit B occurs; 
a new stable orbit K is created and orbit B becomes unstable. Orbits C and C are unstable 
for A30 > and their values of TrM become larger as A30 increases. According to the argument 
in §2.4, we thus expect that the contributions of these orbits to the shell effect decrease with 
increasing A30. On the other hand, orbit A is stable and its TrM value approaches towards 2 as 
A30 increases. This implies that the contribution of orbit A becomes more important. In this way, 
relative magnitude of the amplitude factors between the two families of orbits with the period 
^ 27r/ct;sh and ~ 2Ti/uj±_ changes so that the interference effect between them becomes stronger. 

The above discussion is based on the expression (|2.11|) obtained by the SPA. We should 
note, however, that our classical phase space contains both regular and chaotic regions, i.e., our 
system is a mixed system. As is well known, such a system is abundant in the resonance regions 
where the SPA breaks down, so that the amplitude factors y4„^ should be evaluated by means 
of a more sophisticated method beyond the SPA, e.g., the uniform approximation.il This is an 
interesting future subject, and we expect that the above consideration will remain valid, as long 
as a qualitative feature is concerned, even when nonlinear effects beyond the SPA are taken into 
account. 



6.2. Scaling Properties and Fourier Analysis 



By virtue of the scaling property, Eq. (^4.5| ), the following scaling rules hold 



S,{E) 

m 

ME) 



7 ' 



(6.1) 
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where denotes the degeneracy of orbit 7; (i^ = 1 for a general orbit and for an isolated orbit 
like D or A' in Fig. 0. 

Let us consider the Fourier transform 

P{s) = J dEe''''E-''l^g{E) (6.2) 

of the level density g{E) multiplied by E~^^'^. (The factor E~'^^'^ is attached here to compensate 
for the energy dependence of the amplitude factor A^; see below.) If one insert the exact level 
density g{E) = J2n^{E — En), it becomes 

n 

This quantity can be evaluated with the use of the eigenvalues obtained by a quantum mechanical 
calculation. On the other hand, if we insert the semiclassical level density ( |2.4| ) in ( |6.2| ) and put 
d = 1 appropriate to non-isolated orbits, then we obtain 

p(^i)(s) = P{s) + ^ A„,^(l)e^"^"^/2(5(s - nT^) . (6.4) 

Here P{s) comes from g{E) and has a peak at s = associated with the orbits of zero length. 
On the other hand, the second term on the r.h.s. gives rise to sharp peaks at s = nT^ associated 
with the classical periodic orbits 7 with periods (and their multiple traversals). Note that, 
owing to the scaling property (|6.1|), periods of the primitive orbits are equal to action S'^(l) 
calculated at E = 1. If the trace formula is valid, one expect P("') ~ P^^'"). Thus, we can extract 
information about classical periodic orbits by calculating P^'v^\ Namely, the amplitude factors 
and the Maslov phases of the periodic orbits may be obtained from absolute values and arguments 
of P(^"") (s), respectively. 0' 

Now, let us evaluate the Fourier transform ( |6.3|) . Since the summation is taken over a finite 
number of quantum levels in practice, we introduce the Gaussian cut-off and define a smoothed 
version of it; 

PAsis) = I ds'P(s') fiis' - s)/As) (6.5) 

where f{x) is Gaussian f{x) = exp(— a;^/2). For ( |6.3D and (|6.4|), we obtain 

^ir^(^) = E E-'/'e-'^-fiEjE^^^) , (6.6) 

n 

P^hs) = Pas{s) + E^n,7(l)e^"^-/((s - nT,)/As) , (6.7) 

n,7 
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where i^max = 1/As. 

We calculate the eigenvalues by a matrix diagonalization method with the deformed oscillator 
bases, and use the lower part of the resulting spectrum. Figure [1^ shows the absolute value of 
P^^\s) for A30 = 0.2 ~ 0.4 calculated with -Emax = ^5huJsh{Mo) ■ The loci of the periods of 
classical periodic orbits and their multiple traversals are indicated by arrows in the figures. We 
see nice correspondence between the peaks of P{s) and the periods of classical periodic orbits. 
Almost all peaks can be explained in terms of the classical orbits, indicating that the properties 
of quantum spectrum is characterized mostly by classical periodic orbits. 

Next, let us notice the A30 dependence of P*^'^™)(s). In Fig. we see that the peak at 
s ^ 2n/ujsh decreases, while the peak at s ~ vr/cjsh grows up with increasing A30. Since heights 
of the peaks in P^'^^\s) indicate intensities coming from the corresponding periodic orbits, this 
implies that the contributions from the orbits with the period ^ 27r/u;sh becomes increasingly 
important as A30 increases. The change in relative intensity as a function of A30 between the 



two families of periodic orbit seen in Fig. ^ may be responsible for the enhancement of the 
supershell effect in the reflection-asymmetric SD potential, in accordance with our discussion in 
the preceding subsection. 



6.3. Quantum Signature of Bifurcations 

In order to see how the bifurcations (resonances) of periodic orbits affect the magnitudes of 
the Fourier amplitudes, let us evaluate the heights of the peak as functions of A30 at the periods of 
the classical orbits. As examples, we take the period-tripling and the period-5-upling bifurcations 
of orbit A, which occur at A30 = 0.36 and 0.25, respectively. In Fig. ^ are plotted the calculated 
values of P^'^^\s) as functions of A30 at specific values of s that correspond to the periods of three- 
and five- fold traversals of orbit A. In accordance with the argument given below Eq. (|2.10|) , we 
find that the heights of the peaks indeed exhibit supremes near the bifurcation (resonance) points, 
but little delays are observed for both cases. To account for this delay, it may be necessary to go 
beyond the SPA. 
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6.4. Angular Momentum Decomposition of the Trace Formula 



As our system is axially symmetric, the angular momentum about the symmetry axis is a 
good quantum number. Thus, the level density can be decomposed as g{E) = J2m9{E;m) with 
m denoting the angular momentum quantum number. Let us derive a semiclassical expression 
of g{E;m). Writing the three-dimensional coordinate vector as q = {Q,<f) with Q = {p,z), the 
Green function may be decomposed as 

oo 

G{q\q'-E)= G((Q",^ + 2n7r),(Q',0);E) 

?i=— oo 

= Y. I dMe'('^+2n-)*^G'(Q", Q ] E, M) 

n •' 

oo 

= Yl e'^^GiQ",Q';E,m) , (6.8) 

m=— oo 

where G denotes a Fourier transform of G with respect to ip = if' — (p, and where the Poisson sum 
formula is used in obtaining the last expression. Taking the trace of Eq. ( |6.8| ), one can derive the 
trace formula for g{E] m) in a way similar to 



giE;m) = -2\m j dQpG{Q,Q; E,m) 



g{E; m) + — V] cos{naa{E)/h - (7r/2)/i„«) , (6.9) 



where denotes the action integrals along the two dimensional closed orbits a in the (p, z) plane, 
Tq, the periods, and fia the Maslov phases. Using the symplectic property of the monodromy 



matrix, Wna can be written as 



Wna = det (1 - M^) = 2 - Tr (M^^ , (6.10) 

where Mq, is a 2x2 monodromy matrix. Note that, for symmetric self-retracing orbits, is 
different from appearing in Eq. ( 2.11| ). It is easily seen that, due to the reflection symmetry 



with respect to the 2;-axis, these orbits in the (p, z) plane have periods half of those in the three 
dimensional space. 

Now, for m = 0, a scaling property holds so that we can use the Fourier transformation 
technique. Since the degeneracy of the orbits is zero in the two-dimensional space, we put d = 
in Eq. (|6.6|) . From the above consideration, one expect that the Fourier transform will exhibit 
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peaks, in addition to those corresponding to the periods of closed orbits, also at half of the periods 
of the three-dimensional sjTiimetric self-retracing orbits. The results of calculation is shown in 



Fig. [T^. Again we find a clear correspondence between peaks of the Fourier transform and periods 
of classical orbits. As expected, peaks appear also at half integer times the period of orbit A. 



§7. CONCLUDING REMARKS 

We have found a clear correspondence between the shell structure, i.e., the oscillatory struc- 
ture in the smoothed level density, and the classical periodic orbits for single-particle motions in a 
reflection-asymmetric SD oscillator potential. We have then shown that the supershell effect, i.e., 
an interference effect between two families of the periodic orbits having periods approximately 
2tx /uj_[_ and 27r/u;sh, develops when the reflection-asymmetric deformation increases. This super- 
shell effect is in clear correspondence with the odd-even effect in N^h pointed out in Refs. and 



[ p!0|] . Possible origins of this enhancement phenomena have been pointed out in connection with 
stabilities of the classical periodic orbits. Quantum signature of the period-tripling bifurcation of 
the shortest-period orbit is also pointed out. 

It should be emphasized that our model Hamiltonian system is a mixed system where chaos 
and tori are intermixed; accordingly, period-multipling bifurcations occur, as we have seen, rather 
frequently when the reflection-asymmetric deformation parameter is varied. As is well known, the 
SPA breaks down at the bifurcation points so that we cannot use the Gutzwiller trace formula for 
the aim of calculating the smoothed level density. Instead, by virtue of the scaling property of our 
model Hamiltonian, we have been able to use the Fourier transformation technique to find the 
quantum-classical correspondence. Properties of the Gutzwiller amplitudes have been used only 
as a guide to qualitative discussions. It is an interesting future subject to investigate the problem 
discussed in this paper by using a more sophisticated method, like the uniform approximation, 
which goes beyond the SPA. 
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APPENDIX A: DERIVATION OF EQ. (|2ll| ) 



Semiclassical expression of the Green function in a /-dimensional system is 



Giq", q'; E) = 



27r 



TT 



(Al) 



where S{q",q';E) is the classical action Jp-dq along a trajectory connecting q' and q" with 
energy E. The determinant Ds in the amplitude factor is given by 



D, 



d^S/dq"dq' d^S/dEdq' 
d^S/dq"dE d^S/dE^ 



(A2) 



Let us consider a three-dimensional system (/=3) with axial symmetry, and define an orthogonal 
coordinate q = rj, () for each periodic orbit. We take ^ along the direction of the trajectory 
and f] perpendicular to both ^ and the azimuthal direction (f. Differentiating the Hamilton- Jacobi 
equations 



H{p" = dS/dq", q") = E 
H{p' = -dS/dq', q') = E 



(A3a) 
(A3b) 



with respect to E and using rj = ( = 0, one obtains 



1 _y^ '9^ d'^S 
~ Y dp- dq'/dE ' 



dp[ dEdq[ 



(A4a) 
(A4b) 



where Sxy denotes {d"^ S / dxdy) . If one differentiates ( |A3a|) and ( |A3b|) with ^' and respectively, 
one obtains 







(A5) 



Thus, the determinant flA^ is written as 



See 


Sec 


Seti' 


Sec 


S^'E 











Srj"E 





Sri"ri' 


Sn"c 


S^'E 





S(^"fji 


Sec 



1 1 



Srf'rj' Sri"C 
Sqi'^i S(^ii(^i 



(A6) 
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If we use coordinates t], ip) which are generally not orthogonal, then we obtain 

1 1 



Sijjfi' Sipip 



(AT) 



where J is Jacobian of the coordinate transformation, J' and J" denoting its value at q' and q", 
respectively, and ip = (f" — (f'. Let us evaluate the trace of the Green function in the stationary 
phase approximation. As usual, the action integral along a closed path may be expanded about 
a stationary point q as 



S{q,q;E) = S{q,q;E) + {q-q 



dS dS 

+ 



dq" dq' 



dq"dq" dq"dq' dq'dq' 



(A8) 



The stationary phase condition requires the second term in the r.h.s. to vanish. This is nothing 
but the condition for the trajectory to be periodic, i.e., p" = p'. Taking the axial symmetry into 
account, we can rewrite Eq. ([A8| ) 



as 



1, 



S{q,q-E) = SiE) + -W{Ov' + 



(A9) 



where 



W{0 = det(l - M)Sr,"r^' = (2 - TlM)Sr,"r^' . 



(AlO) 



M being the (2 x 2) monodromy matrix for the periodic orbit (see §2.2). Performing the Gauss- 
Fresnel integral with respect to rj, we finally obtain the following result 

1 



^osc(^) = ^^j^Im^ J dlpd^dr]J^J\Ds\exp 
47r 



{271%^/^ 



E 



2 - TrM^ 



cos(S'^//i — ;U^7r/2) 



where /x^ 



sign(Vr^)/2 and 



I Sr,"r 



1/2 



/ dt 




1/2 


= f dt 




-1/2 


n 






J7 







(All) 



(A12) 
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TABLES 

TABLE I. Properties of the periodic orbits: periods T (in unit of l/a;sh(0)) and traces of the re- 
duced monodronly matrices TrM, cvahiatcd for A30 = 0, 0.2, 0.3 and 0.4. Here, '717' denote the n-fold 
traversal of the primitive orbit 7. For the isolated orbit A' and D, the monodromy matrix M has two 
unit-eigenvalues and the remaining four eigenvalues appear in pairs (e"*,e~"^) and {e°^^ , e~°^^) . These 
pairs are identical (eta = Ob) for orbit D, but they differ from each other for orbit A'. Traces of these 
pairs are given for orbit A'. 



A30 
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0.3 
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FIGURES 

FIG. 1. Results of semiclassical calculation of the oscillating part of the level density for the 
SD oscillator model, a) The contribution from orbits with T = 2t[ /ujz-, b) the contribution from 
orbits with T = 27r/c(j_L, and c) the supershell structure caused by the interference between the 
above two families of orbits. 



FIG. 2. Oscillating parts of the level density for the Hamiltonian ( ^.4| ) with A30 = 0.2, 0.3 and 
0.4, calculated by means of the Strutinsky method with the smoothing width 7 = O.S/icUgh- The 
arrows indicate the minima associated with the even-A^sh closure. 

FIG. 3. The oscillating part of the level density at A30 = 0.4 obtained by the Strutinsky 
method (broken line) is compared with a superposition (solid line) of cos{ET^/h) with the 
period evaluated for the classical orbit A (T ^ 2T[/uj_i_) and that for orbits B, C, C, and 2A 
(T ^ 2-n /ojsh) (see §5 for the properties of these classical orbits). Amplitudes and phases of these 
cosine functions are determined so that the solid line best agrees with the broken line, except that 
the energy dependence of the amplitudes is assumed to fulfill the relation in Eq. (|6.1| ) determined 
by the scaling property of the system under consideration (see §6). 



FIG. 4. Poincare maps in the section {z,Pz) for the Hamiltonian ( |4.1| ) with = and with 
A30 = 0.2 ~ 0.4, defined by = and Pp < 0. 

FIG. 5. Poincare maps in the section {z,Pz) for the Hamiltonian ( |4.1|) with A30 = 0.4 and with 
Pip/E = 0.2 and 0.4, defined hj Pp = and Pp < 0. 

FIG. 6. Poincare maps in the section {p,Pp) for the Hamiltonian (|4.1| ) with = and with 
A30 = 0.2 and 0.4, defined by 2; = and pz > 0. 



FIG. 7. Short periodic orbits for the Hamiltonian (|4.1j ) with A30 = 0.4. Upper part: Planar 
orbits in the plane containing the symmetric axis z. Lower part: A circular orbit in the plane 
perpendicular to the symmetry axis (A') and a three-dimensional orbit (C). Their projections on 
the [x, y) plane and on the {z, y) plane are shown. 
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FIG. 8. Short planar orbits for the Hamiltonian ( |4.1|) with p^p = and A30 = 0.3 ~ 0.4. 



FIG. 9. Traces of the reduced monodromy matrices TrM for the non-isolated periodic orbits 
shown in Fig. |^ (see text for their definitions). 

FIG. 10. Fourier transform P^s^\s) of the level density g{E) defined by Eq. ( |6.6| ), for A30 = 0.2, 
0.3 and 0.4. Gaussian cut off is done with -Emax = h/As = 15/;.ci;sh(A3o)- Arrows indicate periods 
of the classical periodic orbits (see Fig. |^) and of their repetitions. This figure is basically the same 
as Fig. 15 in our previous report; Q-* but, accuracy of the numerical calculation is significantly 
improved so that the peak at s/tt ^ 1 is now clearly seen. This improvement greatly facilitates 
the discussion on the classical-quantum correspondence (see text). 



FIG. 11. Peaks heights in the Fourier transform defined by Eq. (|6.ti|) at the periods of classical 



periodic orbit A and of its multiple traversals (m=l,3,5), plotted as functions of A30. Gaussian cut 
off is done with -Emax = 12/iu;sh(A3o)- Two arrows represent the period-tripling and period-5-upling 
bifurcation points. 

FIG. 12. Fourier transforms of the level density g{E, m) in the m=0 subspace for A3o=0.2, 0.3 
and 0.4. 
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